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' | i Abstract. Quantum mechanics courses focus mostly on its computational aspects. 

This alone does not provide the same depth of understanding as most physicists 
have of classical mechanics. The understanding of classical mechanics is significantly 
, bolstered by the intuitive understanding that one acquires through the playing of 

5j ' games like baseball at an early age. It would be good to have similar games for 

quantum mechanics. However, real games that involve quantum phenomena directly 
are impossible. So, computer simulated games are good alternatives. Here a computer 
game involving three interacting quantum particles is discussed. It is hoped that 
such games played at an early age will provide the intuitive background for a better 
understanding of quantum mechanics at a later age. 
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^ ; PACS numbers: 01.50.-i, 01.50.Wg, 03.65. Ta 
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O . 1. Introduction 

.£h ■ "... and then the wavefunction collapses." What visual images are inspired by such a 

J>Y statement? We can visualize collapsing bridges, buildings and maybe even a souffle. But 
collapsing wavefunctions are a visual mystery for both novices and experts in quantum 
mechanics. A few years of graduate school can teach a physics student the mathematical 
methods as well as the experimental tests of quantum mechanics. However, acquiring an 
intuitive understanding of the subject is more challenging. Classical mechanics is easier 
to understand due to the ready availability of visual images (collapsing bridges, souffles, 
etc.). It is also significant that these images of classical mechanics are observed by 
everyone at an early age making them part of our intuition. Similar early introductions 
to quantum phenomena would be very useful for the learning experience of children. 
They would build a foundation for later, more rigorously mathematical, presentations 
of quantum mechnics. However, real visual images for quantum mechanics are difficult 
to find. So let us look for some computer simulated visual images through a computer 
game based on quantum mechanics. 

Using computers for physics education has become quite mainstream through the 
last decade[TJ 121 EJ E]. However, using physics based computer games is not that 
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Figure 1. The game screen at starting time (left) and at a later time (right). 



common jSJ E]. Personally, I prefer the game approach for several reasons. It can 
provide an intuitive understanding at a very early age without the need for mathematics. 
It is non-threatening and develops physics intuition in a relaxed setting. It is like 
understanding projectile motion while playing baseball. In particular, for quantum 
mechanics, computer games are uniquely useful as real games like baseball shed very 
little light on the subject. 

In the past, I have developed a game based on the quantum mechanical free particle 
( "Quantum Duck Hunt" ) [7J |H] . Here, I present a significantly more complex system - an 
interacting three particle system. The game based on this system is called "Quantum 
Focus" 0. It deals with various subtle aspects of quantum observation and wavefunction 
collapse [TUj. 

This game is not meant to teach quantum mechanics to college students. It 
is meant to develop "quantum intuition" at a much earlier age (maybe elementary 
school). Children playing this kind of games are expected to develop "gut feelings" 
about quantum phenomena just as they usually do for classical phenomena by playing 
baseball or soccer. The approach here is not that of standard accepted padagogy. 
Something different is being tried to see if it works better. I have tried it on a few 
children (my own and their friends!). The results are very encouraging - now they 
want to learn quantum mechanics! If we wait until college to develop quantum intuition 
in youngsters, we might miss the formative years when most intuition is developed. 
However, college students may also benefit from this game by studying the computer 
code and trying to come up with their own variations. 

2. The game 

The game is started by choosing "Start" from the "Action" menu. Three boxes colored 
red, green and blue appear on a black background screen. With time, the color of 
each box begins to smear into neighboring boxes (see figure [TJ. At the same time, the 
brightest spot in each color smear moves away from the other two. As the colors spread, 
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they produce mixtures of the primary colors in various proportions in different boxes 
on the screen. 

The object of the game is to make each color smear as small as possible and at the 
same time bring the three colors as close together as possible. The only tool available 
for achieving this is the click of the mouse at strategic points. Clicking the mouse will 
retract a color completely into a single box - the one that was clicked. But this retraction 
or "collapse" will occur only with a probability proportional to the preexisting intensity 
of that color in that box. So, if there is very little red in a certain box, it is unlikely that 
red can be collapsed into it by any amount of clicking. This is why, sometimes, a click 
of the mouse may produce no effect at all. An interesting sound effect will accompany 
the actual collapse of a color. 

The quantum mechanically minded readers must already have noted that the three 
primary colors represent three particles. The intensity of a primary color is related to 
the probability of finding the corresponding particle at a given place. The mouse clicks 
simulate a particle detector's attempts at detecting a particle. The retraction of a color 
into a single box simulates a wavefunction collapse. In the present model, the three 
particles are tied together by attractive forces. So, it is a bit tricky to see why quantum 
mechanics makes the three color smears move apart. This effect will be discussed later. 

A score is computed in each time step. It depends on how small each color smear 
is and how close the three colors are. So, the goal is to produce a single white box and 
no other colored boxes. But this state of perfection can be seen to be impossible. The 
score displayed (at the top left corner) is the maximum score achieved during the course 
of a game. 

It should be noted that, while the colors spread, nothing is lost. Colors that spread 
off-screen on one edge reappear on the opposite edge. This effect may be used for game 
strategies. 

The game can be played at four levels of difficulty. The scoring formula respects 
the level of difficulty. The features of these levels will be discussed later. 



3. The quantum three particle problem 

This game is based on the quantum dynamics of three interacting distinguishable 
particles. Most quantum problems deal with the solution of the time-independent 
Schrodinger equation. But here we are concerned with the wavefunction collapse and 
subsequent change in the wavefunction. Hence, it is necessary to solve the time- 
dependent Schrodinger equation [TU]. 

Let the positions of the three particles be q i; and their momenta be pi (i = 1, 2, 3.). 
Let the wavefunction of the system be i>({<li},t) and its hamiltonian ff({qj}, {pi}) 
where t is time. Then the time-dependent Schrodinger equation to be solved is fF2\ ITS] 

iM = m (i) 
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Figure 2. Three particles attached by "springs" 



For simplicity I choose harmonic ( "spring" ) potentials to represent the interparticle 
forces (see figure EJ. Also, The springs are assumed to have zero unextended length. 
The resulting hamiltonian is as follows. 

where rrii is the mass and ki the spring constant for the i particle. q c is the position 
of the common center at which the three springs are tied. 

It can be seen that, like most three particle problems, this cannot be separated in 
variables. The dependence on q c complicates the hamiltonian significantly as it is not 
an independent coordinate. q c depends on the particle coordinates due to the following 
zero net force condition. 

3 

- Qc) = 0, (3) 



=i 



which gives 



Hence, 



(4) 



^ = E£- + E^(^-3^J , (5) 




where K = X)i=i %■ 
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Equation ^ can be solved numerically to obtain the time development of the 
wavefunction provided an intial value is specified. At any point in this time development, 
if a particle detector detects the first particle in a small region R, it will be with a 
probability 

P 1 = / d 3 qi / d 3 q 2 / dW({q*},*MH},t)- (6) 

JR J — oo J — oo 

If the particle is actually detected, the wavefunction must collapse to 

MM, t) = AA R (qx - qoMR}, t), (7) 

where A/j(q! — q ) is a sharply peaked function that is nonzero only in the region of 
detection R centered about the position q . The detailed form of this function depends 
on the detector sensitivity in the region R. In the limit R — > 0, it is the Dirac delta 
function: 

lim A R (q fc - q ) = S 3 (q k - q ). (8) 

H — >U 

The constant A is needed to renormalize ip after the collapse. The detection of the other 
two particles can be described similarly. 

After the collapse, ip is replaced by ip c and the time development continued as given 
by equation [T] until the next collapse. 



4. The numerical technique 

The computer screen being 2-dimensional, the above formulation will be reduced to its 
2-dimensional equivalent for the purpose of the game. To use standard finite difference 
methods, the screen space is divided into a matrix of m columns and n rows to produce 
a total ofmxn boxes. For the purpose of a game, we may sacrifice accuracy for speed 
as long as the qualitative aspects of the system are maintained. So, the maximum values 
of m and n are chosen to be 12 and 8. 

To solve the Schrodinger equation, boundary conditions must be specified. There 
are several possible natural choices: 

(i) Perfectly reflecting boundary conditions. 

(ii) Perfectly absorbing boundary conditions. 

(iii) Periodic boundary conditions. 

The perfectly reflecting boundary produces a discontinuity at the boundary that 
interferes with visualization. The perfectly absorbing boundary allows particles to go 
off screen, thus making them useless for visualization. The periodic condition seems 
to be the best for visualization. It identifies the left edge to the right and the bottom 
edge to the top (toroidal topology). Hence, particle current that disappears on one edge 
reappears on the opposite edge. 

The discrete forms for the x and y components of each coordinate q^ may be written 

as 

QlX O'ixAX^ Qiy CljyA?/, (9) 
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where i = 1,2,3, = 0, 1, 2, . . . , m — 1, and = 0, 1, 2, . . . , n — 1. Ax is the mesh 
width in the x direction and Ay is the mesh width in the y direction. If the mesh 
width in time is At, then equation ^ produces the following recursive formula for the 
computation of if). 



In general, the above numerical algorithm for the solution of the Schrodinger equation 
is known to be unstable [TT]. However, we can use it in the present case because 
wavefunction collapses are expected to preempt any instability. Besides, as noted before 
we are not looking for high accuracy. 

The wavefunction if), at one instant of time, is a function of all coordinates q^. So, 
its discretized form must depend on all ai x and (%. Thus, for numerical computation, ip 
is represented by an array of 6 dimensions (one for each aj X . and a^). In the notation of 
the C language it would be: , 0[ai a: ][ai y ][a2x][a2?/][a3a:][a3i;]- For compactness of notation 
I can write this as: if) [a] [b] [c] [d] [e] [f] or ifj a ,b,c,d,e,f- Then the finite difference form of 
the operation by the hamiltonian H is found from equation 03 using equation E3 and the 
following finite difference forms of the operators. 



Here the most common finite difference form for second derivatives is used. Using 
equations El El and[TT]in equation [TU1 the wavefunction for successive time steps can be 
computed. The numerical method chosen here does not maintain normalization of i/j. 
Hence, after each time step computation, ip must be normalized (TIJ. 

Also after each time step computation, the screen image must be updated to provide 
an animated visual effect. The RGB coloring scheme on the computer screen provides 
a natural way of representing the three particle probabilities. The red, green and blue 
color intensities in a box are made proportional to the probabilities of finding each 
of the three particles in that box. The discretized version of equation |H1 provides the 
probabilities to be used. They follows. 




(10) 





c,d,e,f 
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P 3 = Yl ^aAcAeJ^AcAeJ- (12) 
a,b,c,d 

So, the C++ code used to find the amount of red color (particle 1) in a box is as follows. 

int CQFocusDoc : :Red(int a, int b) 
{ 

int c , d , e , f ; 
float red = 0; 

f or (c=0 ; c<xpts ; C++) 

for (d=0 ; d<ypts ; d++) 

f or (e=0 ; e<xpts ; e++) 

for(f=0;f<ypts;f++) 

red += norm(Psi [a] [b] [c] [d] [e] [f ] ) ; 

return (255*sqrt (sqrt (red) ) ) ; //Square root used to enhance 

//color for better visibility. 

} 

The amounts of the other two colors are computed similarly. 

When the mouse button is clicked in a box, one of the three particles is picked 
randomly for collapse and then the decision to actually collapse it is made based on 
the probability given by equation El The following C++ code fragment makes these 
probabilistic decisions. 

partnum = MyRandom(3) ; //Generates an integer random number between and 2. 
intensity = MyRandom(256) ; 

intensity = (intensity*intensity*intensity) / (256*256) ; 

// The above redifinition improves game by requiring less mouse clicking. 

i = point .x/bwidth; // x pixel position divided by box width, 
j = point .y/bheight; // y pixel position divided by box height. 

swit ch (partnum) 
{ 

case 0: 

if(intensity < pDoc->Red(i, j)) // Function Red(i,j) defined above. 

{ pDoc->Collapse (partnum, CPoint(i, j)) ; setcollapse = true;} 
break; 
case 1 : 

if (intensity < pDoc->Green(i , j) ) 

{ pDoc->Collapse (partnum, CPoint(i, j)) ; setcollapse = true;} 
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break; 
case 2: 

if (intensity < pDoc->Blue(i, j)) 

{ pDoc->Collapse (partnum, CPoint(i, j)) ; setcollapse = true;} 
break; 

} 

The wavefunction after the collapse is given by equation [7| The function A^(q^ — q ) 
in its discrete form is chosen as the discrete form of the Dirac delta function: 



where the integer values a^x, clo x , o-ky and a 0y are defined as in equation^ So, the C++ 
code for wavefunction collapse is as follows. 

void CQFocusDoc :: Collapse (int partnum, CPoint point) 
{ 

int a, b, c, d, e, f; 

switch (partnum) 

{ 

case 0: 

f or(a=0;a<xpts; a++) 
for (b=0 ; b<ypts ; b++) 
f or (c=0; c<xpts ; C++) 
for(d=0;d<ypts;d++) 
f or (e=0 ; e<xpts ; e++) 
for(f=0;f<ypts;f++) 
{ 

if (a!=point .x I I b!=point.y) 
Psi[a] [b] [c] [d] [e] [f] =0; 

} 

break; 
case 1 : 

f or(a=0;a<xpts; a++) 
for (b=0 ; b<ypts ; b++) 
f or (c=0; c<xpts ; C++) 
for (d=0 ; d<ypts ; d++) 
f or (e=0 ; e<xpts ; e++) 
for(f=0;f<ypts;f++) 
{ 

if (c ! =point .x I I d!=point.y) 
Psi[a] [b] [c] [d] [e] [f] =0; 

} 




(13) 



break; 
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case 2: 

f or(a=0;a<xpts; a++) 
for (b=0 ; b<ypts ; b++) 
f or (c=0; c<xpts ; C++) 
for(d=0;d<ypts;d++) 
f or (e=0 ; e<xpts ; e++) 
for(f=0;f<ypts;f++) 
{ 

if (e ! =point . x I I f!=point.y) 
Psi[a] [b] [c] [d] [e] [f] =0; 

} 

break; 

} 

Normalize () ; 

} 

Here partnum gives the randomly picked particle number and point identifies the box 
clicked. 

5. The levels of difficulty 

The game can be played at four different levels of difficulty. The difficulty level is 
increased by increasing the values of the spring constants This increases the rate at 
which the positions of maximum probability move apart. The reason will be seen in the 
next section. 

The difficulty level is increased also by reducing the total number of boxes. This 
increases the speed of computation and hence increases the rate of spreading of the 
wavef unction. 

The score allows for higher difficulty levels. 

6. Some results 

The primary purpose of this game is to provide repeated and consistent visual effects 
that mimic quantum wavefunction dynamics. As expected, the wavefunction collapse 
leaves the undetected particles unaffected. Also as expected, the probability profile 
of each particle spreads with time. The resulting mix of the primary colors produces 
some rather unusual color effects that may interest the artists amongst us. What is 
not-so-obvious is as follows. If we start with one particle in a collapsed state (with no 
velocity), with time its probability peak moves away from those of the other particles! 
As the potential function used here is attractive, this is somewhat surprising. However, 
closer scrutiny can explain this phenomenon. 

Consider the standard one-particle harmonic oscillator. Higher energy eigenstates 
have probability peaks farther away from the origin. This means that particles that 
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start off with higher momenta are likely to have their probability peaks farther away. 

For the present case, a particle collapsed to its position eigenstate has high probabilities 

for large momenta and hence, large energy. This makes its probability peak move away 

from the other particles. 
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